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Revisiting Zeno's paradox: Flying arrows for atom-diatom reactions 
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The possibility to observe quantum Zeno and anti-Zeno scenarios for atom-diatom reactive colli- 
sions is investigated for two diferent processes (F-f HD and H+O2) by means of time-dependent wave 
packet propagations. A novel approach is proposed in which the survival probabilities investigated 
are those obtained when the initial state is redefined after time intervals r at which measurements 
are performed on the system. The comparison with the actual probability for the unperturbed 
system reveals the existence of a regime in which the decay from the initial state appears to be 
inhibited (the quantum Zeno effect) or accelerated (the anti-Zeno effect). In contrast with preced- 
ing interpretations, given that the time-evolving wave packet is not affected at any time {^^ flying 
arrow"), the present procedure does not require to invoke to counterintuitive hindered evolutions 
{''^stopping arrow") to explain the results. On the contrary the consecutive update of the reference 
state of the system solves the apparent modern version of Zeno's paradox. 

PACS numbers: Valid PACS appear here 



In the extreme limit of one of his paradoxes, Zeno pro- 
posed many centuries ago the abolition of motion. The 
Greek philosopher defied the common sense when he con- 
cluded that a flying arrow was at rest at each point of 
time during its motion to the target. In order to solve 
the apparent controversy, one is then forced to invoke the 
differential calculus by defining the velocity of the arrow, 
not by its instantaneous position in space, but via the 
corresponding time derivative of the spatial coordinates. 
The spirit of this postulate has been adopted nowadays 
in the investigation of the evolution of a quantum unsta- 
ble system when frequent measurements are performed 
on it. According to the quantum Zeno effect (QZE) and 
anti-Zeno effect (AZE), one can either inhibit or accel- 
erate, respectively, the decay of a specific quantum state 
by means of observations on the system. These phenom- 
ena have their origin in the deviations from the expected 
exponential decay of the corresponding survival proba- 
bility at the short-time regime [T]. Various experimental 
confirmations on a broad list of subfields in Physics have 
been reported In particular, the occurrence of both 
effects has been observed in an ensemble of cold sodium 
atoms trapped in an accelerated standing wave of light 
[3]. The authors found that the decay on the number of 
remaining trapped atoms depends on the frequency of the 
measurements, and it could be suppresed or enhanced in 
comparison with the unperturbed system. 

From a theorical point of view, Kofman and Kurizki 
[4] established that the modifications introduced to the 
decaying dynamics of the observed system are due to 
the energy spread induced by the measurements and the 
spectral density of the states coupled to the unstable ini- 
tial state of the system. Facchi et al. [5] have investi- 
gated the transition between the QZE and AZE regimes 
in terms of the interval separating successive measure- 
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ments. The intersection between the survival probabiHty 
and the exponential decay of the state of the correspond- 
ing system according to its "natural" rate, 70, defines the 
value of the time interval, r*, for consecutive measure- 
ments which actually separate both regimes. 

The Zeno effect has also been investigated for chemi- 
cal reactions. In particular, Prezhdo [6 found evidences 
of quantum AZE in the photodissociation of the van der 
Waals (vdW) cluster HF-Na. The process can be accel- 
erated by an order of magnitude by controlling the loss 
of electronic decohcrcncc induced by interactions between 
the system and the enviroment. The choice of the specific 
solvent could then be used to monitor the fragmentation 
process. In a similar spirit, we have recently studied the 
competition between the vibrational and electronic pre- 
dissociation mechanisms for the NeBr2 vdW cluster [7]. 
The branching ratio between both fragmentation chan- 
nels was found to be significantly modified by slowing 
down the fast vibrational predissociation pathway with 
successive measurements on the system. 

In this letter we address the occurrence of the Zeno 
effect in the context of the reaction dynamics of atom- 
diatom collisions studied by means of a time-dependent 
wave packet (WP) approach. The time propagation can 
be written as: 

|^(t))=e-^^*/''Uo), (1) 

where describes the value of the WP at i = 0. The sur- 
vival probability or decay function from the initial state 
is simply given by: 

p{t) ^ {cp, I m p . (2) 

With the measurement, the unitarity of the process, 
as described in Eq. ([I]), is broken down and the quan- 
tum state describing the system "collapses" into one of 
the pointer states, (/)o, associated with the measuring de- 
vice. In particular, as a result of an ideal selective and 



2 



complete measurement [5] to investigate the decay of the 
system at consecutive r intervals, the state of the system 
at i = Nt becomes: 

I ^Nr) = [(</>0 I ^r)f I M, (3) 

where the prefactor is a consequence of the projector 
P =1 </)o)(0o I- The corresponding survival probability 
can then be expressed as: 

PNrit) =1 (<^0 I <t>r) P^l {<t>^ I m) ? ■ (4) 

The comparison of the corresponding survival proba- 
bilities expressed in Eq. Q for different intervals t and 
the unperturbed probability P(t) (see Eq. (l2j) reveals 
the possible occurrence of both QZE and AZE |5j . Those 
Pnt (t) probabilities with a slower decay in time with re- 
spect to P{t) are usually interpreted as delayed evolu- 
tions of the system when we perform measurements to 
analyse its state, in close analogy to the premises stated 
in Zenos's paradox {"stopping arrovJ^ approach). The 
evolution of the WP is governed by a constantly in- 
terrupted sequence of time operators which act within 
the time intervals between the measurements performed 
a.t t — nr. 

However it is possible to conceive a different rationale 
to explain the behaviour of the survival probabilities. 
Inspired by one of the key remarks made in Ref. [3], 
regarding the definition of a new initial state after the 
measurement, we can devise an alternative approach in 
which the state of the system at each time t = t inter- 
val, is employed as a reference to analyse the evolution of 
the time-propagating WP. Thus, the survival probability 
after Nt measurements is written as follows: 

PNrit) =1 i^iN-Dr I ^A^r) P| {<pNr \ </>W> |' • (5) 

The prefactor in Eq. ([5| can be initially understood as 
the constant which ensures that the survival probabil- 
ity reinitiates at the value taken a.t t = {N — 1)t dur- 
ing the time interval [{N — 1)t,Nt]. In this scheme, the 
WP (j){t) does not alter its time evolution {'^flying arrow" 
approach) which is governed by the expression given in 
Eq. ([!]) and therefore, it is possible to express (I)nt and 
<P(n-i)t in Eq. ([s]) as a function of according to Eq. 
(n]). Thus, the survival probability is expressed as: 

PNr{t) -I (^/-O I <i^r) n {<l>Nr \ Ht)) l' • (6) 

This latter approach, in which the initial WP employed 
as a reference to trace the time evolution of the system 
is redefined, but the evolving WP, (j){t), as in the case of 
the flying arrow in Zeno's paradox, is not perturbed, has 
been highlighted by Fischer et al. [Sj in order to detect 
the QZE and AZE. In their experiment, authors of Ref. 3 
adopted the remaining number of atoms in the trap after 



the measurement as the initial condition. The presence 
of the prefactor in Eq. ([6]), originated by the projection 
performed at each measurement, thus plays an analogue 
role to the ansatz invoked in the experiment of Ref. ^ to 
reset the time to t = 0. 

In this letter we discuss this new approach for two re- 
actions: F-I-HD and H+O2. The time dependent WP ap- 
proach employed here has been described before [8, and 
it is based on the description of the total Hamiltonian, 
H = T + V in terms of mass-scaled Jacobi coordinates. 
The propagation is carried out by means of the symmet- 
ric split-operator method [S] and the application of an 
absorbing potential [8] for the outgoing fluxes at limiting 
regions of both reactants and products arrangements to 
avoid numerical reflections. For the corresponding V of 
the F-I-HD reaction, we used the three dimensional SW 
[To] potential energy surface (PES) while the DMBEIV 
[TP PES was employed for the H-I-O2 reactions. The ini- 
tial WP for the F-I-HD case is located at the transition 
state, around the region of the PES at which a reactive 
resonance at a collision energy of about 0.52 kcal/mol 
was found [12 . For the H-I-O2 reaction, the initial WP 
started its propagation at the asymptotic region of the 
reactants channel. 

The corresponding survival probability P{t) as defined 
in Eq. ([2| is shown in Figure [l] (black solid line) for both 
processes. The curves exhibit a pronounced decay with 
a minimum for t ~ 0.5 fs, for the F-I-HD reaction and 
t ^ 1 is for H-I-O2, due to a fast dephasing between the 
initial WP and 4>{t). The natural exponential decay be- 
yond the short-time regime, Ze"'*'"* [5], is obtained from 
the tail of the secondary maximum features developed at 
t = 0.9 and 1.7 fs for the F-I-HD and H+O2 reactions 
respectively. Both, the Ze^"'°'^ curve and the e~''°* prob- 
ability which describes the exponential decay of the sys- 
tem from t = with this natural decay rate, are shown, 
with dashed lines in Figure [l] In the discussion of Ref. 
[5] regarding the existence of a r* value which separates 
the QZE and AZE regimes, it was found that for Z > 1, 
the curves e"'''"' and P{t) may not cross. Interestingly, 
each of the two reactions under analysis in this work, cor- 
respond to a specific case: For the F-I-HD, we find that 
Z > 1, whereas for H-I-O2, Z < 1. 

The e"^"* curve and the unperturbed P{t) probabil- 
ity are found to intersect at r* ^ 0.06 fs for the F-I-HD 
reaction and 0.32 for the H-I-O2 reaction. We have cal- 
culated the survival probability according to Eq. ^ for 
two different values of time intervals between consecutive 
measurements such that t < t* (0.013 fs and 0.026 fs for 
the F-KHD case and 0.066 fs and 0.132 fs for the H+O2). 
These Pt(^) decay probabilies, shown in Fig. [T] extend 
noticeably longer in time than the original P\t): they 
can be then interpreted as a manifestation of the QZE. 

Alternatively, the same procedure is applied to mea- 
surements spaced in time a value t > r*. In particular, 
we have considered r = 0.132 fs for F+HD and 0.527 fs 
for H-I-O2. The crossing between the corresponding Prit) 
survival probabilities and the unperturbed P{t) curves 



3 




FIG. 1: Survival probability P{t) for the H+O2 (left panel) and F+HD reaction (right panel) when no redefinition of the initial 
WP is performed (solid black curve), and when this is done at different t — t values. The curve with the natural decay rate 
70 of the system observed beyond the short-time regime Ze and the decaying probability that the system would present 
if no deviations from the exponential behaviour were observed e~^°* are included with dashed line. The crossing between the 
unperturbed P{t) curve and e"^"* at r* is shown with a vertical dotted line. 



for each reaction finally manifests much faster decays to 
zero when we operate on the system, an evidence of the 
AZE. 

It is worth mentioning that the survival probabilities 
obtained after the operational procedure described in 
Eqs. ([5]) and ^ exhibit a profile in which the recur- 
rence observed for the unperturbed P(t) probabilities is 
completely lost. Those PNrit) probabilities, on the other 
hand, exhibit an oscillating structure with values at each 
t = TIT time interval which follow perfectly, as discussed 
in Ref. 5i the exponentially decay governed by the 7(t) 
rate. 

Whereas the possible differences observed on the sur- 
vival probabilities obtained as in Eq. ^ with respect 
to the unperturbed probability P{t) could certainly be 
explained in terms of variations of the evolution of the 
system from its initial state (at least in the subspace de- 
fined by the measurement [2]), the interpretation for the 
results obtained via the new approach proposed in this 
letter involves more subtle considerations. Under the as- 
sumption of the projected state after the measurement 
as the new starting point for the propagation, the actual 
evolution of the system, here represented by the time 
dependent WP, is not affected at all: the WP remains 



unaltered and no differences in the observed dynamics is 
thus expected after the measurements on the system. In 
particular, the observed QZE could be understood, not 
as a hindered evolution, but as a consequence of the up- 
date of the reference used to measure the overall change 
of the state of the system. 

On the other hand, the occurrence of the AZE seems to 
be subject to some geometrical considerations, as stated 
by Facchi et al. [S]: the definition of r*, the limiting 
value between the QZE and AZE of the time intervals 
for the measurement, is based on the expected exponen- 
tial decay of the survival probability. In this sense, the 
probabilities shown in Fig. [T] exhibit such an exponen- 
tial behaviour once a secondary maximum is developed, 
in a similar way as the examples discussed in Ref. [5] 
The T* value is then required to be sufficiently large to 
ensure that the P{t) function is low enough at the mo- 
ment of the measurement; thus the ulterior evolution will 
yield probabilities which do not exceed the limit imposed 
by the exponentially decaying behaviour followed beyond 
the short-time regime. 

Our findings are consistent with a short-time dynam- 
ics. First of all, the time scale of Fig. [TJis in line with this 
specific regime, with most of the interesting features oc- 
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curring well below 2 fs for both the H+O2 and F+HD re- 
actions here investigated. Furthermore, one might won- 
der if the survival probabilities as expressed in Eqs. Q 
and ^ are, in essence, formally equivalent, despite of 
the fact that the corresponding ip{t) WPs follow different 
time evolutions. This seems to be the case if we could ap- 
ply the definition of the ^^vt function given in Eq. ^ to 
expand the PNrit) probability of the flying arrow" ap- 
proach (see Eq. The redefined initial state at each 
t = T would be then a rescaling of the state of the system 
at t — 0, 00, with a pref actor which reduces the corre- 
sponding amplitude. The above mentioned short-time 
regimes seem to be the restriction for the strict validity 
of such approaches. 

Since our analysis is enterely based on the projection of 
the initial state and their corresponding time-updated 
expressions (j)r, one might wonder about the possible uni- 
versality of the occurrence of the QZE and AZE scenarios 
here described. In our calculations an ample set of dif- 
ferent locations and average incident momenta, po, for 
the initial WP were found to lead to similar results, but, 
for the F-I-HD reaction, the most favourable cases were 
those in which the WP was located near to the transition 
state, the choice discussed in this paper. 

The dependence of the occurrence of both QZE and 
AZE on the parameters of the initial WP can be better 
analyzed in the one-dimensional case. For an initial state 
described by a Gaussian WP in free space, V{x) — 0, 
the survival probability in the short time limit can be 
expressed as [13]: 



Pit) 



-2{po/p,f(t^/8tl) 



(7) 



where Ps = h/2ao is the spreading momentum and tc = 
2to(To /fi^ is a characteristic time with ctq being the initial 
spreading of the WP. If the translational and spreading 
motion of the WP are such that pa/ps ^ '^tc/t, then Eq. 
([7]) is written as: 



P{t) = l- 



Po 



(8) 



which displays the characteristic quadratic behaviour 
with time of Zeno's scenario [2]. For any other relation 
between pq and Ps, the survival probability decays too 
fast to observe the QZE. In the above Eq. ([s]), if the 
translation momentum of the WP is too small to be ne- 
glected in comparison with the spreading motion, po « 0, 
we find that the r* time which separates the QZE and 
AZE is given by r* = 2^/2 1^. 

According to this, the QZE will be observed whenever 
the measurements are performed before the WP has en- 
tered the time-scales of linear spreading. If, on the con- 
trary, the WP has some non-zero initial momentum, we 
find that: 



(9) 



Our present study revisits Zeno's paradox for time- 
propagating WPs in the context of atom-diatom reac- 
tions. In order to understand the effect of frequent mea- 
surements on the survival probabilities, we propose a 
novel approach: the redefinition of the initial state used 
as a reference to calculate the decay of the state of the 
system. Thus, apparently counterintuitive explanations 
regarding the evolution of the time evolving WPs are not 
required. A possible extention to the present study could 
include resonances and quasibound states, examples of 
time-decaying features within the context of molecular 
processes treated here. 
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